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ABSTRACT 

Spitzer mid- infrared images of the dusty warped disk in the galaxy Centaurus A show a parallelogram-shaped struc- 
ture. We successfully model the observed mid-infrared morphology by integrating the light from an emitting, thin, 
and warped disk, similar to that inferred from previous kinematic studies. The models with the best match to the 
morphology lack dust emission within the inner 0.1 to 0.8 kpc, suggesting that energetic processes near the nucleus 
have disturbed the inner molecular disk, creating a gap in the molecular gas distribution. 

Subject headings: galaxies: structure - galaxies: ISM - galaxies: individual (NGC 5128) - galaxies: 
peculiar 



1. INTRODUCTION 

Centaurus A (NGC 5128) is the nearest of all the giant 
radio galaxies. Because of the disk of gas and dust in its 
central regions, Centaurus A is suspected to be the prod- 
uct of a merger of a small gas rich spiral galaxy wit h 
a larger elliptical galaxy ijBaade fc Minkowskil 119541) . 
Numerical simulat ions of such mergers pr edict large 
shell-like features ijHernauist fc Quinnlll988T) that have 
been observed in Centaurus A over a la rge range of 
radii (|Malin et all 119831: iPeng et all 12002(1 . Some con- 
tain atomic and molecular gas, implying that they 
were formed from material stripped from a galaxy rel- 
atively re cently, fewer than a few galactic rotation pe- 
riods ago jSchiminovich et al.lll994HCharmandaris et al.l 
2000). Tidal features associated with this debris have 
also been identified, supporting the relatively short (few 
times 10 8 years) estimated timescale since the merger 
took place l|Peng et alJl2002fh 

In its central regions, NGC 5128 exhibits a well rec- 
ognized, optically-dark band of absorption across its 
nucleus. This dusty disk w as first mode led as a 
transi e nt warped st r ucture bv iTubbsl 11980J). j Bland 
llT98l. iBland et all jpsl - iQuillen et all l)1992fl . and 
iNicholson et alJ l|1992|) found that the kinematics of 
the ionized and molecular gas are well modeled by 
a warped disk composed of a series of inclined con- 
nected rings undergoing circular motion (as also ex- 
plored for other galaxies wit h peculiar morphology by 
iSteiman-Cameron et al.lll992fl . The model explored by 
Quillen, Graha m fc Frogel (1993) modified the kine- 
matic model bv IQuillen et al.1 (fl992) to fit the morphol- 
ogy of the absorptive, dusty disk seen in near-infrared 
images and proposed a timescale of about 200 million 
years since the core of an infalling spiral galaxy reached 
and merged with the elliptical galaxy nucleus. An ini- 
tially flat disk, misaligned with the galaxy principal sym- 
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metry axis, becomes increasingly corrugated as a func- 
tion of time. The short timescale estimated since the 
merger in NGC 5128 is approximately consistent with 
the timescale suggested by the presence of tidal debris 
and by the shell-like features containing ato mic hydrogen 
llSchiminovich et alJll994HPeng et al]l2002T) and molecu- 
lar gas ( Charmand aris et al.l2000I) . An alternative model 
ac counting for th e warped disk is the polar ring model 
bv ISparka l|1996j) . consistent with the polar orbit of the 
disk implied by the galaxy isophotes in the outer galaxy 
lIMafin et al.lli983t lHavTles et a,UlT983l IPeng et al]l2002ft 
but requiring a longer timescale to account for the twist 
of the warp. 

More recent observations of the central region include 
submillimcter imaging with the Submillimeter Com- 
mon User Bolometer Array (SCUBA) and mid-IR imag- 
ing with the Infrared Camera (ISOCA M) on the In- 
frared Space Obse rvatory (ISO) satellite (Mirabel 1999; 
iLeeuw et al.ll20 02). At these wavelengths, the dusty disk 
is seen in emission rather than absorption. At 100-200 pc 
from the nucleus the galaxy contains a molecular cir- 
cumnucle ar disk that has b e en stu died in molecular line 
emission ( [Israel et al.lll990t 119911) and resolved in Pa a 
emission (Sch reier et al.lll99fiHMarconi et alJ l2001). For 
a recent summary of the wealth of observational stud- 
ies carried out on this peculiar a nd active nea rby galaxy, 
see the comprehe nsive review bv llsT^lp9^ . Based on 
the discussion bv llsraej l(1998t) . we adopt a distance to 
Centaurus A of 3.4 Mpc. At this distance 1' on the sky 
corresponds to ~ 1 kpc. 

In this manuscript we focus on the geometry of the 
dusty disk in Centaurus A as seen from Spitzer Infrared 
Array Camera (IRAC) images. These images resolve the 
structure of the disk more clearly than previous obser- 
vations and also show the disk out to larger radii. They 
provide us with an opportunity to study the geometry 
of the warped disk in much higher detail than previously 
possible. Observations are described in § 2. Our geomet- 
ric model is described in § 3. A discussion and summary 
follow in § 4. 

2. OBSERVATIONS 
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Figures and present images of NGC 5128 taken on 
2004 February 10 in the 3.6, 4.5, 5.8 and 8.0 /jm broad- 
band filters (cha nnels 1-4) of the S yitzer Infrared Array 
Camera (IRAC: iFazio et all 120041) . In each filter, fixed 
cluster observing mode was used for these observations 
to produce a map of 5 x 6 points, at which five dithered 
exposures were taken. The exposure time per frame was 
12 s. Additional shorter exposure, 0.4 s frames, were also 
taken to correct for possible saturation in the longer ex- 
posure frames. The coverage of the map at each position 
on the sky varies between three and ten frames, with 
an average coverage of six frames, corresponding to an 
exposure time of 72 s. 

Before mosaicing, the basic calibrated data (BCD) 
frames were corrected for artifacts using the IRAC arti- 
fact mitigation code (excluding the pulldown correction) 
available from the Spitzer Science Center contributed 
software pages 1 . The final mosaiced images were pro- 
duced from the dithered frames by applying the MOPEX 
software 2 to the corrected BCD frames. 

Here we focus on the central 5' x 5' region in the final 
mosaics. The mosaiced images provide a plate scale with 
1"2 pixels; FWHMs of the point spread functions are 
1.7, 1.7, 1.9, 2'.'0 in channels 1-4, respectively. The rms 
noise levels in these images are 0.025, 0.024, 0.070, and 
0.060 MJysr -1 in channels 1-4, respectively. The above 
sensitivities agree with the predictions of the SENS-PET 
sensitivity estimator in longer wavelengths channels. The 
3.6 /im channel contains emission from the stellar com- 
ponent of the galaxy that extends over a large portion 
of the field of view. This makes the 3.6 /im image about 
half as sensitive as the SENS-PET prediction for the dust 
emission. 

The IRAC images pr esented here ar e deeper than pre- 
vious ISOCAM images l|Mirabell l|1999j) ). and have higher 
a ngular resolut i on than previous s ubmillimeter images 
( ljMirabelllT999t ILeeuw et al1l2002|) ). In the inner few 
arcminutes, the IRAC images exhibit a parallelogram 
shape in emission (see Fig. 0). This shape, previously 
seen at lower angular resolution, was int erpreted as an 
S-shape, possibly a ssociated with shocks (Mirabel 1999; 
ILeeuw et al.ll200l . A parallelogram morphology has 
been seen previously in other galaxies. For example, dust 
absorption features in the SO galaxy NGC 4753 exhibit 
a parallelogram shape, and have been modeled with a 
warped twisted disk by Steima n-Cameron et al.l (1992). 
When an optically thin warped disk is seen in emission, 
the edges of folds in the disk correspond to regions of 
higher surface brightness. Mul tiple folds are seen along 
the line of sight (e.g.. IBlandl Il986l: IBland et all 119871: 
IQuillen et aTlll993t INicholson et alJl!992|) implying that 
some parts of the disk are seen through other parts of the 
disk. In this case the morphology would not be nearly 
symmetric across the origin (r — » — r), as observed, un- 
less the outer disk was nearly optically thin in the mid- 
infrared. We infer that the disk probably has a low nor- 
mal optical depth in the mid-infrared (if observed face- 
on), though the bright edges of the folds and individual 
clumps in the disk may not be optically thin. The par- 
allelogram is present in all four IRAC bands, though it 
is most prominent compared to the diffuse emission from 

1 http: / /ssc. spitzer. caltech.edu/archanaly/contributed/browse. html 

2 http:/ /ssc. spitzer. caltech.edu/postbcd/ 



starlight in the longer wavelength 5.8 and 8fim bands. 

In visible or near-infrared bands, folds can correspond 
to regions where the absorption from dust obscures back- 
ground starlight. In this case folds that are closer to 
the observer absorb more background starlight from the 
galaxy. Thus some of the features seen in emission in 
the IRAC images resemble and coincide with absorption 
bands previously seen in near-infrared images. In Fig- 
ure [3] we show a color map made from 2 Micron All Sky 
Surve y (2MASS) images from the 2MASS large galaxy 
atlas ijJarrett et all 2003). The south-eastern edge of the 
parallelogram in the IRAC images lies in the same loca- 
tion as the absorption band ~ 10" to the south east of 
the nucleus prominent in the n ear-infrared images . This 
is consistent with the study bv ILeeuw et all ((2002), who 
compared submillimeter images to near-infrared images. 
The edge of emission ~ 1' to the north of the nucleus 
(at a surface brightness below that of the parallelogram) 
in the 8.0 /zm IRAC image corresponds to the top of the 
dust lane prominent in optical images of the galaxy. This 
corresponds to the extinction feature in the near-infrared 
color map to the north of the nucleus. From comparing 
the IRAC images to the 2MASS images we infer that the 
southern side of the parallelogram is closer to the ob- 
server than the northern side. Because the dust on this 
side lies in front of the plane perpendicular to the line 
of sight containing the galaxy nucleus, it absorbs more 
background starlight and so causes a deeper band of ex- 
tinction in the near-infrared images. The northern side 
of the parallelogram lies on the opposite side and so is not 
seen in the near-infrared images. Likewise the oval edge 
of emission ~ 1' to the north of the galaxy nucleus (see 
Fig. |2J) is nearer the observer than that on the opposite 
edge ~ 1' to the south of the nucleus, and only the north- 
ern side causes an absorption feature in the near-infrared 
images and color maps (see Fig- • 

3. MODELING OF THE WARPED DISK 

At many observed wavelengths, the morphology of 
Centaurus A is well r eproduced by geometric mod- 
els of a warp ed di sk (IBlandl 119861 IBland et all 119871: 
I Quillen et aTTll99a INicholson et alJl!992t IQuillen et all 
119931 ISparkdll996|) . We describe and extend such mod- 
eling here. 

A warped disk can be described as a series of tilted 
rings, each with a different radius r. We assume that 
the gas and dust are evenly distributed on each ring and 
each ring is smoothly connected to those at larger and 
sm aller radii. We follow the no tation and framework used 
by IQuillen et all <|1992l Il993j) to describe the geometry 
of a warped disk. Each ring is described by two angles, 
a precession angle a(r) and an inclination angle uj(r). 
These angles are given with respect to an assumed prin- 
cipal axis of the underlying elliptical galaxy. This axis 
requires two angles to describe, Xj corresponding to the 
position angle counter clockwise from north of the axis 
on the sky, and an inclination angle, that describes 
the tilt of this gal actic axis with resp ect to the line of 
sight (see Fig. 6 bv lQuillen et al.lH993j) . We assume that 
the galaxy is axisymmetric and not triaxial so a third 
angle is not required to describe its orientation. We also 
assume that the galaxy shape is fixed, and not tumbling. 

Our description for the ring projection angles can be 
related to those of previous works. While IQuillen et alJ 
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l)1992tll993|) described the angles of the warped disk with 
respect to the princ ipal axis of the underlying el l iptica l 
galaxy, iBland et alJ <|1987f> and INicholson etafl (|1992f) 
matched the Ha velocity field with a tilted ring model 
describing the orientation of the rings with respect to 
the line of sight. This model fit the velocity field by 
adjusting the ring inclination as a function of radius, i(r), 
and position angle counter clockwise from north, p(r). 

To produce a model image of the mid-infrared emis- 
sion, we must consider all emitting and absorbing regions 
along the line of sight at each position on the sky. This 
is much simpler when the system is optically thin. In 
this case, we sum all emission along the line of sight at 
each position on the sky. The near symmetry of the disk 
suggests that we can use an optically thin approxima- 
tion to model this disk. While individual clouds could 
contain optically thick regions, the bulk of the emitted 
mid-infrared light is likely to reach the observer. When 
the disk is optically thin, brighter areas correspond to 
regions that appear folded from the perspective of the 
viewer. Here we neglect emission from a spherical com- 
ponent associated with the stars and only consider emis- 
sion from a thin but warped disk. 

Our numerical procedure begins by randomly sampling 
x, y positions in the plane perpendicular to the principal 
axis of the elliptical galaxy. At each position we compute 
a z-coordinate based on a smooth (spline) function for 
the disk precession and inclination angles, a(r) and u>(r). 
The points that specify these spline functions are listed in 
Table 1. The coordinates for each point are then rotated 
into the viewer's frame. Points along each line of sight 
are summed to produce a model image. For the functions 
a(r) a nd w(r), we began with those from lOtiillen et all 
(1993) and slightly varied the angles at different radii 
to achieve a better match to the observed morphology. 
Matching was done by visual comparison to the IRAC 
images. No minimizing fit to the image data values was 
done. Our procedure is adequate to understand the pro- 
jection affects associated with an emitting warped disk. 
Future and more intensive modeling would be required 
to do a multi-dimensional and multi- wavelength fit. 

To construct a model, we must consider the thickness 
of the disk and its brightness distribution. Because the 
disk is not infinitely thin we randomly chose slight off- 
sets from the z-coordinate computed from the precession 
and inclination angles. The size of the offset depends 
on a disk aspect ratio k(r) = h/r, where h(r) is the disk 
thickness as a function of radius. We assumed for the disk 
aspect ratio a power law form, fc(r) = kso(r/5Q")P* . The 
intensity of the emission contributed at each point was 
then computed assuming that the disk volume emissiv- 
ity integrated through the disk vertically is a power law 
function of radius, e(r) oc r~^ c . We also allowed an in- 
ner gap in the radial surface brightness profile, r gap . The 
contribution from the active nucleus and the ~ 100 pc cir- 
cumnuclear disk are not taken into account in our model 
so by a gap we mean a deficit within r gap and that of the 
circumnuclear disk at ~ 6". 

A model intensity image computed as described above 
is shown in Figure 0] along with the 8.0 fim IRAC image 
of the galaxy. The numerical parameters of this model 
are summarized in Table ^ The angle precession and 
inclination angular functions, a(r) and i(r), are shown 
in Figure as a solid and dashed lines, respectively. The 



precession angle used previously bv lQuillen et alJ l)1993fl 
from matching the near-infrared morphology is shown 
as a dotted line in this figure and is similar to our best 
matching model. The Spitzer data allow us to better 
study the outer parts of the disk than possible with the 
old 12 CO(2-l) spectra and near-infrared images. Thus 
there are some differences in the precession angle a(r) 
between our old model and our newer one, as shown in 
Figure El To better match the morphology in the outer 
parts of the disk we allowed the disk to twist further (to 
higher a higher value of a) before decreasing at larger 
radius. 

Our model predicts that the disk alternates between 
having the southern side and northern side nearest, with 
folds both above the nucleus. The near-infrared images 
show a strong dust absorption feature to the southwest of 
the nucleus corresponding to an inner fold, and a weaker 
featu re north of the nuc leus corresponding to an outer 
fold ijQuillen et al.lfl99^ . In Figure we show the rings 
comprising our warp model, projected on the sky. The 
nearer semi-circle of each ring is shown in red, whereas 
the more distant semi-circle is shown in blue. The disk 
at r ~ 60" is nearest the observer on the south-east side. 
This region corresponds to a fold in the disk that is seen 
in the near-infrared color map, 10" to the south-east of 
the nucleus (see Fig- - At r > 100" the northern side of 
the disk is again nearest the observer. This corresponds 
to the northern dust lane seen in the the optical images 
and the near-infrared color map about 45" north of the 
nucleus. 

3.1. Comparison to previous studies on the geometry of 
the warped disk 

In Figure 01, b, we show the inclination and position 
angles as a function of radius with respect to the viewer; 
fu nctions that can be dir ectly compared to those found 
bv lNicholson et all l|1992|) from their kin ematic fits to the 
Ha ve locity field (see their Fig. 9a,b). INicholson et alJ 
(1992) adopted a distance of 3 Mpc to the galaxy, so their 
linear scale is 15% different than ours; we have corrected 
for this difference by rescaling the given distances. We 
see from Figure [7| that the position angles of their tilted 
ring fit are quite similar to those predicted by our model. 

The warp shape originally d esigned to match the ve- 
locity field in molecular gas l|Quillen et alJ Il992fl re- 
sembled that which matched the Ha velocity field by 
INicholson et al. (1992). Both models included an in- 
clined disk that tilted so that gas rings at different radii 
alternated between being retrograde and prograde with 
respect to the observer. This was also a characteris- 
tic of the high inclination warp model for NGC 4753 
by Steiman-Cameron et al.l l)1992|) . The tilted ring fit 
by INicholson et al.l l|1992fl to the velocity field did not 
specify which side of th e each r ing was clo s est th e 
viewer. However iBlandl II1986D. IBland et all |1987ft . 
INicholson et all (|1992|) . and iQuillen et all l)1993|) used 
the visible and near-infrared images to break this degen- 
eracy. The disk is tilted so that it is nearer the viewer on 
the southern side for intermediate radii; accounting for 
the southern edge of the dust lane seen in optical images, 
and nearer the viewer on the northern side at largest radii 
to account for the northern dust lane seen in the optical 
and near-infrared images. These flips in orientation can 
be seen from Figure [TJo. They occur where the inclina- 
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tion with respect to the viewer crosses 90° at which point 
the disk is edge-on. 

The model by iQuillen et al.l l)1993|) had a somewhat 
lower inclination for the galaxy principal axis, 1? = 65°, 
instead of 75°, used here. Their model also had a higher 
disk inclination with respect to the galaxy principal axis 
at smaller radii. There is redundancy in the model be- 
tween the inclination of the disk with respect to the 
galaxy axis, w(r), and the tilt of the galaxy principal 
axis, The maximum and minimum ring inclinations 
correspond to valu es set by i? ± w (see discussion by 
IQuillen et al"]ll992fl . To exhibit the parallelo gram shape, 
the disk inclination with respect to the viewer must go 
above and below 90° causing folds in the disk. This re- 
stricts the models to a range of values $ + ui ~ 100°. 

Here, we find a somewhat better match (see Fig. [HJ to 
the mid-infrared images with a model that has decreasing 
inclination (with respect to our estimated galaxy princi- 
pal axis) at larger radii. However this decrease may be 
spurious. To compute the model at sufficient resolution 
to compare to the inner region, we cannot well sample 
points over a large region. By accounting for every ob- 
served feature we may have achieved a better morphology 
match with a model that has exaggerated variations in 
the ring angles. In other words, the angular variations 
may have been compressed into a smaller radii than they 
should be. 

To be more certain of the orientation of the outer 
disk, better kinematic constraints on this outer disk 
are needed. By specifying a choice for the principal 
axis of the galaxy, we have chosen to describe the ge- 
ometry of the warped disk with respect to a pa rticu- 
lar axis. This choice helped IQuillen et alJ l) 19931) com- 
pare the shape of the disk to the predictions of sim- 
ple merger and precessing ring models. However, it 
is not necessarily significant that the disk inclination 
as measured with respect to our assumed galaxy axis 
varies slowly with radius. An additional complexity 
not considered with our choice of projection angles is 
that the galaxy may be triaxial, or significantly vary 
in shap e and orientation with rad ius (e.g., as consid- 
ered by lArnaboldi k, Sparkd" Il994|) . The galaxy could 
be in the process of dynamic relaxation following the 
merger. Analytical models or one-dimensional integra- 
tions fail to capture the complexity of more detailed nu- 
merical s imulations, particular for near equal ma ss merg- 
ers (e.g.. lMihoslll999t iMihos fc Hernauistlll99fil) . Imag- 
ing and kinematic studies of the outer galaxy suggest 
that these additional degrees of freedom are important 
llMalin et all 11983 ISchiminovich et alJ 119941 iHui et all 
ll995HPeng et al.ll2002t) ~ 

3.2. In context with the dynamical warp models 

The model precession angle a(r) reaches a maximum 
at r ~ 100", decreasing at both larger and smaller radii 
(see Fig. |SJ). This implies that the corrugated disk is 
twisted in one direction for r < 100" and in the op- 
posite direction for r > 100". Since the near-infrared 
morphology depends on which side of each ring is closer 
to the observer, and our model is similar to that used 
to match the near-infrared images (also se e discussion 
bv lBland etaD 119871 IRicholson et atlll992tt . this pecu- 
liar change in the handedness of the twist is likely to be 
real. When the handedness of the twist is the same at all 



radii, a nested se t of parallelogram shap es can be seen 
(e.g., as in M84; IQuillen fc Bowerl ll999). However in 
CentaurusA, an outer oval is seen in the outer disk (see 
Fig. |S]and[3l that is slanted in the opposite direction as 
the inner parallelogram. This difference in the direction 
of the slant is a feature of the sign change in the slope of 
the precession angle. 

The precession rate of a gas ring inclined with respect 
to the underlying galaxy undergoing circular motion is 
proportional to the ellipticity of the galactic gravitational 
potential and the angular rotation rate of the ring (e.g., 
lTubbdll980l ISoarkdll98^ . If the self gravity of the ring 
is important t hen it too can affect the precession rat e 
llSparkel 11984 lArnaboldi fc Sparkd IT9941 ISparkel Il99^ . 
The direction of precession depends on the orientation 
of the ring and galaxy (whether a polar ring or not) and 
whether the galaxy is prolate or oblate. Most models 
assume that following a merger the gas and dust are dis- 
tributed in a plane that is misaligned with the principal 
axes of the galaxy. To be consistent with the direction of 
the twist in the outer parts of Centaurus A (decreasing 
a(r) at r > 100"), the galaxy can either be prolate and 
the r ing located near the equatorial plane ijQuillen et alJ 
11993ft or the galaxy can be oblate and the ring would be 
nearly polar l)SparkdlT 996). 

The reversal in the sign of the slope (da/dr) is unlikely 
to be caused solely by the shape of the rotation curve. 
The predicted rotation curve of Centaurus A increases at 
s maller radius; all the w ay to a radius of 10" (see Fig. 13 
bv lMarconi et al. 1120 05) . Consequently the angular rota- 
tion rate should increase as the radius decreases nearly all 
th e way to the g a laxy c enter The rotation curves used 
by IQuillen et~aLl (1993) and ISparkd l|1996|) were nearly 
solid-body at small radii, and so underestimated the an- 
gular precession rate at small radii (within ~ 1.5' of the 
nucleus) . 

The reversal in the sign of the slope (da/dr) could 
be due to a drop in the galaxy eccentr icity (exploited 
by both IQuillen et alJll993USparkdll996D that would re- 
duce the precession rate at small r adii. To account for 
the change in slope of a at r ~ 100", Qui llen et alJ l|1993fl 
assumed a sharp cutoff in the ellipticity of the gravita- 
tional potential of the galaxy as a function of radius. 
For r < 80", the potential ellipticity was much reduced 
in the model by Quille rTet al.l l)1993|) compared to that 
outside it. The galaxy isophotes are best viewed in the 
short wavelength 3.6/^m IRAC image where the stellar 
light most contributes to the flux and the extinction and 
emission from dust is minimized compared to that at 
shorter wavelengths in the near-infrared. Galaxy isopho- 
tal contours for this image are shown in Figure [5J 

We first examine the isophotes at 3.6/im to search 
for evid e nce o f the large scale stellar bar proposed by 
iMirabel |[l999). If such a bar were responsible for the 
parallelogram shape in the mid-infrared images, the bar 
would be viewed at intermediate inclinations (not edge- 
on), and so should be evident in the isophote shapes. 
However, the isophotes do not exhibit a region of flat 
surface brightness or a change in ellipticity or position 
angle over a short range in radius. These features would 
be expected at the end of a stellar bar. We conclude 
that there is no evidence for a large scale, few kpc sized, 
stellar bar at the heart of CentaurusA. 

The galaxy isophotal ellipticity does decrease in the 
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inner regions. By fitting contours to the 3.6 /im image 
outside the emitting disk, we measured a galaxy isopho- 
tal ellipticity for the stellar distribution (e = 1 — b/a, 
where b and a are the semi-minor and semi-major axes) of 
e ~ 0.1 at r ~ 4', e ~ 0.05 at r ~ 2', and e ~ at r < 1'. 
Because the gravitational potential is a convolution of 
the density profile with al/r function, the gravitational 
potential contours are smoother than the associated den- 
sity distribution. However, the isophotes are likely to be 
rounder and smoother than the actual density distribu- 
tion. While the drop in background ga l axy ellipticity 
exploited bv lQuillen et all l|1993h : lSparkel l|1996ft in their 
dynamical models is real, it may not be sharp enough 
to account for the abrupt drop in the precession rate in- 
ferred in the central region of the galaxy from the change 
in the slope of the precession angle a at r ~ 100". 

We now co nsider the role of the mass in the disk. 
iSparkd 1)1996(1 showed that the self gravity of the disk 
could increase the precessi on rate and va ry the precession 
axis in the central region. iSnarkel l)1996t) only considered 
the mass in the mo lecular and atomic gas components. 
Quill erTet alJ l|1993f) noted that there were extensions in 
the K-band isophotes that were not reproduced by the 
purely absorptive disk model. The K-band isophotes are 
extended at a radius of about 60" in the 15 mag/arcsec 2 
contour (als o see the dereddened contour map shown in 
Figure 9 bv lMarconi et al] l20051. We can compare the 
mass in a possible disk gas and stellar component to that 
in the underlying galaxy. With a circular rotational ve- 
locity of 250 km s -1 , a total mass of ~ 10 10 M Q is en- 
closed within a radius of 1 kpc. The mass in molecular 
gas in the same region is a few times 10 8 M© (Phillips 
et al. 1987, corrected for the difference in the assumed 
distance). The level of 15magarcsec~ 2 in K-band cor- 
responds to a surface density of 4500 Mq pc~ 2 assuming 
a ma ss to light ratio at K-band of 0.5 (e.g., Sil ge et aLl 
2005). The actual surface density would be ~ 10 times 
lower than this, taking into account the mean disk axis 
ratio on the sky. This surface density can be compared 
to the estimated surface density of molecular gas, or a 
few hundred -M©pc~ 2 . This suggests that at least an 
equal mass exists in stars in the disk as in gas. 

Comparing the total mass in gas and stars in the disk 
to that in the underlying spherical component, we esti- 
mate that a few percent of the total mass within 1 kpc 
lies in the disk. The mass in th e disk is likely to be a few 
times larger than that used bv ISnarkel l)1996l) to account 
for the disk geometry. We support the finding by Sparkc 
(1996) that the self-gravity in the disk is important, and 
so should significantly affect the disk precession rate. Fu- 
ture modifications of the prolate model should take this 
into account, as the self-gravity in the disk could change 
the direction of precession in this region. If the outer 
disk is prolate, then a reversal in the twist could be due 
to the self-gravity of the disk in the inner region. If the 
outer gal axy is oblate, then the ring is polar and the 
model by Sparke (1996) would account for the reversal 
in the twist of the disk. Both dynamical models could 
be updated to include better estimates of the mass in 
the disk, the galaxy isophotes at 2.2 and 3.6 /im, and an 
improved rotation curve based on the light distribution. 

3.3. A gap in the dusty disk between a radius of 6" and 

50" 



The model that is most similar to the IRAC images 
(shown in Fig. 0} contains a gap in the dust distribution 
with outer radius r ~ 50". We compare the model shown 
in Figure 01 to a similar one that lacks the deficit in the 
dust distribution. This model is shown in Figurc lTUb . and 
does not match the observed morphology as well as that 
containing a gap in the dust distribution. The smooth 
continuation of the precession angle into the inner region 
results in an edge-on disk with respect to the viewer at 
some point within r — 60". This causes the sharp bright 
linear feature at r < 60" seen in Figure HUh that is not 
exhibited by the IRAC images. 

If the disk has a lower inclination with respect to the 
viewer, the observed surface brightness is reduced. We 
consider the possibility that the disk precession angle dif- 
fers in the inner region from that expected from a smooth 
continuation of our model. To test this possibility, we 
computed a model that has a flattened (non-twisted) disk 
at r < 60". This model is shown in Figure HHb. The 
sharp, bright, inner feature along the east- west direction 
seen in Figure lTUk is not evident in Figure lTUb . However, 
the morphology of this model also does not display the 
characteristics of our preferred model shown in Figure 0] 
that better resembles the IRAC images. In particular, 
the triangular features to the south-east and north-west 
of the nucleus that are see n in the near-infrar ed color 
maps on the south-east side l)Quillen et alJll993T) are not 
as good a match to those of the IRAC images. To re- 
move the bright inner feature of Figure 110b . the preces- 
sion angle must remain above a > 270°. However, the 
curved triangle edge of emission in the parallelogram is 
not present if the precession angle does not steeply drop 
between a radius of 100" and 50". We have tried mod- 
els with both increasing and decreasing precession angles 
near the nucleus, finding no improved match to the ob- 
served morphology. Our best match is the model with an 
inner gap in the dust distribution. We explored models 
with a lower but constant surface density within the es- 
timated gap radius of ~ 50", and can exclude those with 
a surface density in the gap that is above 1/5 of that at 
the outer gap edge. 

Previous studies have discussed the possibility of a 
deficit in the gas distribution in the same region as we 
find a deficit in the dust distribution. At 100-200 pc from 
the nucleus there is a circumnuclear molecul ar disk that 
has been st udied in molecular line emission l|Israel et alJ 
ll990Lll99lll . A deficit in the ionized gas distribution out- 
side t he circumnuclear disk was seen in the Pa a kinemat- 
ics bv Marconi et all l|2001|) . In Pa a, the circumnuclear 
disk at a radius of r ~ 6" is seen, and so is emission at 
significantly lower velocities (hence inferred larger radii) 
and significantly higher velocities (within the sphere of 
influence of the massive black hole) . However Pa a emis- 
sion is lacking at intermediate radii and velocities. This 
lack of emission would be expected if there were a deficit 
of gaseous material between the radii of ~ 6" (set by 
the estimated oute r radius of the circumnuc lear disk, 
ISchreier et al l liol iMaTcTim et al .1 1200(1 l200l and 50" 
(estimated from our model). 

The rotation curve previously fit to the CO and Ha 
kinematics rose lin early (solid body) w i thin a radius of 
1' of the nucleus l)Quillen et al.l Il992t iNicholson et alJ 
1992). Because an edge-on gas ring appears linear on a 
position-velocity diagram, a gas disk with an inner hole 



() 



can mimic or be confused with a gas disk extending all 
the way to the nucleus in a galaxy wi th a linearly rising 
rotation curve. iNicholson et alJ l)1992f) showed that there 
was a discrepancy between their measured rotation curve 
and that expected from a r 1 / 4 or deVaucoulers law. They 
listed a possible hole in the HII region distribution as a 
possible cause for this discrepancy. A linearly rising rota- 
tion curve within a radius of 1' of the nucleus is not con- 
sistent with the K-band surface brightness profile (see the 
rotat ion curve predicted in Figure 13 bv lMarconi et all 
I2005D . The apparent region of solid body rotation need 
not be account ed for with a galactic bar, as proposed by 
lMirabe]llfl999h . and can be better explained by a gap in 
the gas and dust distribution. 

4. SUMMARY AND DISCUSSION 

By integrating the light through an emitting, opti- 
cally thin, dusty, and warped disk, we have successfully 
matched the morphology of CentaurusA seen in mid- 
infrared IRAC images. We confirm previous proposals 
that the disk morphology is well explained by a warped 
disk iBlandllT98?llBland ' et alJH987t IQuillen et al.lll998t 
lNicholso^i*eT^dT( l^^L^ u^et^!T 20021) rather than a 
barred one l)Mirabellll999[K The disk is nearly edge-on 
with respect to the viewer, but tilts so that folds ap- 
pear above and below the galaxy equator. There is a 
fold south of the nucleus at a radius r ~ 60" responsible 
for high extinction seen in near-infrared images, and a 
fold north of the nucleus at larger radii responsible for 
the northern edge of the dust lane seen in optical and 
images. Extinction features seen in the near-infrared ex- 
tinction map correspond to folds in the disk that are 
located nearer the observer and so absorb more back- 
ground stellar light from the galaxy. In the mid infrared, 
however, folds on both the near and opposite side of the 
galaxy correspond to bright emission features. The in- 
ner folds account for the parallelogram shape, while the 
outer folds account for the northern edge of the dust lane 
seen in optical images and a fainter oval of emission seen 
outside the parallelogram in the IRAC images. 

The disk geometry we use to match the mid-infrared 
morphology i s sim ilar to that found previously by 
IQuillen et alJ l)1993|) . and is also similar to that re quired 
to fit the CO and Ha velocity field iBland et all Il987t 



lillen et aTll!992t INicholson et al.lll992l). The geomet- 
ric wa rp models by Guillen ct al. ( 1993); Nicholso n et al.l 
( 1992) have been predictive; they provide good matches 
to the mid-infrared morphology. Some differences in the 
precession angle exist in the model at r > 100" com- 
pared to previous work. Previous CO, Ha spectra and 
near-infrared imaging were not sensitive enough to pro- 
vide tight kinematic constraints on the outer disk. The 
Spitzer data have allowed us to extend the model past 
r ~ 100" compared to previous models and see closer 
in to the nucleus where the extinction is high at shorter 
wavelengths. Better constraints on the disk geometry 
could be achieved in the future by fitting observations at 
more than one wavelength. For example, modern high 
resolution kinematic observations that are fit along with 
the Spitzer data using the same model would allow much 
stronger constraints on the disk geometry and gas and 
dust distribution than the slight modification to previ- 
ous geometric models that we have used here. 
The warp disk model suggests that there is a gap in 



the dusty disk at 6" < r < 5 0". A gap exists in sa me re- 
gion in the gas distribution. iMarconi et alJ l|200lD saw a 
deficit of Pa a emission near the nucleus at intermediate 
velocities. INicholson et alJ Jl992) suggested that the dis- 
crepancy between the measured linearly rising rotation 
curve within a radius of 60" might be explained by a hole 
in the ionized gas distribution. It is not easy to deter- 
mine if there is a gap in the dust distribution since the 
infrared surface brightness depends on the inclination of 
the disk and the disk is highly corrugated. If the disk 
twists to lower inclinations (closer to face-on) at small 
radii, it would have a lower surface brightness near the 
nucleus. We have explored models with smoothly vary- 
ing radial surface brightness distributions, and smoothly 
varying precession and inclination angles. We find that 
only models with a deficit of dust interior to ~ 50" re- 
semble the mid-infrared images. We exclude models that 
have more than ~ 1/5 the surface density within r = 50" 
as at that radius. We conclude that there is a gap in the 
gas and dust distribution between 0.1 and 0.8 kpc from 
the nucleus. The inner radius of the gap we have taken 
from studies of the ci rcumnuclear disk l)Marconi et alJ 
l2001lllsrael et al.ll991|) . and the outer radius is estimated 
from our model. 

It is interesting that only the region between the cir- 
cumnuclear disk (100-200 pc) and ~ 0.8 kpc has been 
depleted of gas and dust. While we find no evidence 
for a large scale stellar bar in CentaurusA, a dense gas 
disk could have exhibited dynamical inst abilities deposit- 
ing g as into the circumnuclear disk (e.g.. lShlosman et alJ 
1989). Energetic star formation or activity associ- 
ated with the black hole can depl ete and evacuate the 
central region of a galaxy (e.g., ISpringel et alJ 120051 
iVeilleux et al .1120051 . though it is not clear how a cir- 
cumnuclear disk would be protected from or reformed 
following this activity. Nearby galaxies exhibit evacu- 
ated central regions or circumnuclear rings. For exam- 
ple, the Circinus galaxy has a ~ 500 pc radius molecu- 
lar ring l)Curran et al.lll99^ and contains a Seyfert nu- 
cleus. M82 contains a 1 kpc radius circumnuclear molecu- 
lar ring; however, a previous epo ch of star formation has 
occur red within this ring (e.g., iForster Schreiber et alJ 
2003) . Future observational studies may differentiate be- 
tween the role of star formation, dynamical instabilities 
and nuclear activity in disrupting the gas and dusty disk 
in Centaurus A. 

The best matching geometric warp model requires a 
change in the slope of the precession angle at a radius 
of about r ~ 100". Tw o prev ious models account for 
this twist. IQuillen et alJ <) 19931) used a model in which 
the galaxy was prolate in its outer region and the galaxy 
ellipticity abruptly dropped to zero within r ~ 80". How- 
ever, this abrupt drop is not consistent with the isophotc 
shapes at 3.6 /im. A prolate model (with the advantage of 
a relatively short timescale) might account for the disk 
geometry if modified to includ e the self-grav ity of the 
disk. The polar ring model by iSnarkel l)1996l) naturally 
accounts for the reversal, but requires a longer timescale 
to operate. These dynamical models could be updated 
to use a more accurate mass distribution and rotation 
curve. Improvements to these models may lead to better 
understanding of the galaxy merger that created Centau- 
rus A's peculiar morphology, as well as the merger's role 
in feeding the active galactic nucleus. 
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Right Ascension 

Fig. 1. — a) The central region of the image at 3.6 fim taken by the IRAC camera on board the Spitzcr Space Telescope. The emission 
is shown on a log scale. At 3.6 and 4.5 jim starlight peaking near the galaxy center is seen in addition to emission from dust. At longer 
wavelengths emission from dust is primarily seen. 



TABLE 1 

Model parameters 



Description 


Parameter 


Value 




Galaxy principal axis 


PA on sky 


X 


20° 


Tilt 


i? 


75° 


Dust intensity 


Dust emissivity index 


Pe 


4.0 


Aspect ratio 


fcso 


0.1 


Aspect ratio index 


Pk 


0.9 


Radius of inner hole 


1"gap 


50" 



Note. — The precession angle for this model, a(r), is a spline function that interpolates between 10 points in the format (radius in 
arcseconds, a in degrees): (0.2, 180), (36.0, 210 ), (57, 270), (80, 345), (89, 375), (100, 400), (135, 380), (155,345), (240,280), (400,240). 
The inclination angle ui(r) is a spline function that interpolates between the following points (0.2, 25), (100, 20), (260, 8), (400, 8). 
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Right Ascension 



Fig. 2. — The same region as Fig, HI but a color composite of the IRAC images at 3.6, 5.8 and 8.0/jm. This figure is shown on a log 
scale. The black and white version of this figure only shows the emission at at 8.0 /im. 
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Fig. 3. — Color map showing the log of the H-band divided by the J-band 2MASS images by Jarrctt et al. (2003). Lighter shadding 
corresponds to regions of heavier extinction. Folds in the disk both nearer and more distant the observer are seen in the mid-infrared 
images. However folds nearer the observer absorb more background starlight and so are more prominent in the near-infrared images. 
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Fig. 4.— 



On the left is the IRAC 8.0/xm image. On the right is the model for the warped disk with parameters given in Table 1. 
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Fig. 5. — This image shows rings for the warped disk model shown in the previous figure. The semi-circle of each ring nearest the 
observer is shown in red, that more distant shown in blue. For the black-white version of this figure: The semi-circle of each ring nearest 
the observer is shown in black, that more distant shown in gray. At near-infrared and visible wavelengths dust nearer the observer absorbs 
more background galactic starlight than dust more distant from the observer. The inner fold corresponds to the extinction feature south- 
east of the nucleus seen in the J/H-band color map (see FigEJ)). The outer fold seen here corresponds to the extinction feature north of 
the nucleus that is the northern edge of the dustlane prominent in visible images of the galaxy. Regions where rings are close together 
correspond to regions of higher surface brightness in the mid-infrared images. Though the nearer side can cause deeper extinction features 
in the near-infrared images, both sides of the disk cause emission features in the mid-infrared images. 
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Fig. 6. — The solid line shows the precession angle, a(r) as a function of radius used to make the model shown in Figure^ The dashed 
line shows th e adop ted inclination angle. For comparison we show as a dotted line the precession angle used for the previous model by 
Omllen^^ah ^^23) (based on the near-infrar ed imaging). Th ese models should be approximately consistent with the CO velocity field by 
Quillen et all 119921) and H a velocity field bv lNicholson et all 119921) . 



13 




250 



r (arcsec) 




250 



(arcsec) 



Fig. 7. — a) The position angle of each ring on the sky (anti-clockwise f r om N orth) is shown. These angles are compared with those 
estimated from tilted ring fits from the Ha velocity field bv lNicholson et alj 119921) that are shown as diamonds. These points were taken 
from their Figure 9a. b) We show the inclination angle of each ring with respect to the viewer for the model shown in the previous figures. 
The inclination of an edge-on disk is shown as a dotted line at 90° . Rings with inclinations above 90° are oriented with their nearest points 
to the north-east of the nucleus, and those with inclinations below 90 ° have nearest poin t s to the south-west of t he nucleus. These angles 
arc also compared with those from the Ho velocity field (diamonds). Bl and et all fl.987); Nicholso n et alj 119921) infer red that at r < 60" 
the front face of the disk is oriented to the south-west of the nucleus, and vice versa for r > 60" (also see lSparkelll996ti . 
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Fig. 8. — Similar to Figure l4l except shown on a larger scale. 
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Fig. 9. — The IRAC Band 1 (3.6/jm) image shown with isophotal contours. The lowest contour is at 0.63MJysr~ . The contours are 
separated by 0.5 magnitudes or a factor of 1.585. The galaxy is more highly elliptical at larger radii. 




-100 100 -100 100 

foresee) foresee) 



Fig. 10. — a) This model is similar to that shown in Figure HI except there is no deficit in the dust distribution within r gap . b) This 
model is similar to that shown in Figure HI except the precession angle a(r) remains above 270° for r < 50". The IRAC images are best 
matched with a model that contains a deficit in the dust distribution within r ~ 50". 



